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Abstract. In these lecture notes, we present the equations presently used in stellar 
interior models in order to compute the effects of axial rotation. We discuss the hy- 
potheses made. We suggest that the effects of rotation might play a key role at low 
metallicity. 

1 Physics of rotation 

Axial rotation modifies the hydrostatic equilibrium configuration by adding a 
centrifugal acceleration term in the hydrostatic equation, induces many insta- 
bilities driving the transport of angular momentum and of chemical species in 
radiative zones and changes the mass loss rates. In the present lecture notes we 
shall consider the case of models without magnetic fields. 

1.1 Hydrostatic effects of rotation 

The equations of stellar structure. In a rotating star, the equations of stel- 
lar structure need to be modified |23j . The usual spherical coordinates must 
be replaced by new coordinates characterizing the equipotentials. The classi- 
cal method applies when the effective gravity can be derived from a potential 
= <P—^f2 2 r 2 sin 2 9, i.e. when the problem is conservative. There, is the grav- 
itational potential, Q the angular velocity, r the radius at the colatitude 9. If 
the rotation law is shellular (i.e. such that Q is constant on isobaric surfaces see 
below), the problem is non-conservative. Most existing models of rotating stars 
apply, rather inconsistently, the classical scheme by [23]. However, as shown by 
[44] . the equations of stellar structure can still be written consistently, in term 
of a coordinate referring to the mass inside the isobaric surface^) . Thus, the 
problem of the stellar structure of a differentially rotating star in a shellular 
rotation state can be kept one-dimensional. 

The Roche model. In all the derivations, we shall use the Roche model, i.e. we 
approximate the gravitational potential by GM-jr where M- is the mass inside 
the isobaric surface with a mean radius r. The radius r which labels each isobaric 

1 For shellular rotation, the shape of the isobaric surfaces are given by the same expres- 
sion as the one giving the shape of the equipotentials in conservative cases provided 
some changes of variables are performed. 
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surface is defined by r = (V- / (4 / 3tt)) 1 where V- is the volume (deformed by 
rotation) inside the isobaric surface considered. Apart from the case of extreme 
rotational velocities, the parameter r is close to the average radius of an isobar, 
which is the radius at P2(cosz9) = 0, namely for z9 = 54.7 degrees. 

In the frame of the Roche model, the shape of a meridian at the surface of 
the star (which is an isobaric surface) is given by couples of R and 9 values 
satisfying the following equation: 

R 2 Rp' V ' 

where R is the radius at colatitude 9, Q the angular velocity, M the mass inside 
the surface and R p , the polar radius. Thus the shape of the surface (as well 
as of any isobaric surface inside the star) is determined by three parameters 
M, Q and R p . The first two M and Q are independent variables. The third 
one is derived from the first two and the equations of stellar structure. Setting 

x = (tjt-)~ 1/3 R, one can write Eq.Q](see [23]) 

- + ^x 2 sm 2 9 = — . (2) 

With this change of variable, the shape of an equipotential is uniquely deter- 
mined by only one parameter x p . 
Setting 

(3) 

Hp 

where R e is the equatorial radius, one easily obtains from Eq. [T]that 
The above equation relates the inverse of the oblateness / to R p . 



The von Zeipel theorem and its consequences. The von Zeipel theorem 
[60] expresses that the radiative flux F at some colatitude § in a rotating star is 
proportional to the local effective gravity </ e ff . [23] has generalized this theorem 
to the case of shellular rotation and the expression of the flux F for a star with 
angular velocity fl on the isobaric stellar surface is 
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There, p m is the internal average density, \ = 4acT 3 /(3Kp) and \t is the partial 
derivative with respect to T. The quantity is defined by — |, i.e. the 
ratio of the horizontal density fluctuation to the average density on the isobar 
[64 . One has the thermodynamic coefficients 5 = —{d\sxp/d\D.T)p tli , Ht is the 
temperature scale height. The term £($), which expresses the deviations of the 
von Zeipel theorem due to the baroclinicity of the star, is generally very small, 
(cf. [29]). 

Let us emphasize that the flux is proportional to g e ff and not to <7tot- This 
results from the fact that the equation of hydrostatic equilibrium is ^f- = —g e ff . 
The effect of radiation pressure is already counted in the expression of P, which 
is the total pressure. We may call M* the effective mass, i.e. the mass reduced 
by the centrifugal force. This is the complete form of the von Zeipel theorem in 
a differentially rotating star with shellular rotation, whether or not one is close 
to the Eddington limit. 

This theorem has numerous consequences. Some of them are discussed below. 

The position of the star in the HR diagram 

A fast rotating star has stronger radiative fluxes at the pole than at the 
equator. Therefore the position of such a star in the HR diagram will depend of 
the angle between the line of sight and the rotational axis (inclination angle). 
If for instance that angle is 90 degrees, a great part of the light will come from 
the equatorial belt characterized by lower radiative fluxes and cooler effective 
temperatures, while when the star is observed pole-on most of the light will come 
from the hot polar region characterized by stronger fluxes and higher effective 
temperatures. Thus the perceived luminosity and effective temperature (and also 
effective gravity) of a star depend on the inclination angle. This has to be kept 
in mind when comparisons are made with observed quantities. Computations 
of the effect of the inclination angle on the emergent luminosity, colors and 
spectrum have been performed by [35j . The effect of the inclination angle on the 
determination of the effective gravity is discussed in [3T] . In general, a theoretical 
evolutionary track is given in term of total luminosity and of an average effective 
temperature defined by T c 4 ff = L/(aS(f2)), where a is Stefan's constant and 
S(f2) the total actual stellar surface. The total luminosity (corresponding to 
the integrated flux over the surface) does not depend on the angle of view, 
but cannot be directly compared to the "observed" luminosity deduced from the 
apparent luminosity coming from the hemisphere directed toward us. Let us note 
however that for surface velocities inferior to about 70% of the critical velocity 
these effects remain quite modest. As a numerical example, the ratio (T c ir(pole)- 
T c ff (equator) )/T e s (equator) becomes superior to 10% only for uj > 0.7. At break- 
up, the effective temperature of the polar region is about a factor two higher 
than that of the equatorial one. 

The Eddington luminosity 

Let us express the total gravity at some colatitude t9, taking into account the 
radiative acceleration (cf. f2§\ ) 
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thus one has 



9tot = 9eff + 9rad = 9eff H (9) 

The rotation effects appear both in g e ff and in F. We may also consider the 
local limiting flux. The condition g to t = allows us to define a limiting flux, 

FiimO?) - -^ySeff W . (10) 

From that we may define the ratio .Tfo(i9) of the actual flux (see Eq. |5J) F{&) to 
the limiting local flux in a rotating star, 



r «W = p — T^T = 7 7^~T ■ ( n ) 



2-KGp 

As a matter of fact, Fq^) is the local Eddington ratio and 

4 ^GAf(l-3^-) 



'Edd 



[1 + CW] 



(12) 



is the local Eddington luminosity. For a certain angular velocity on the isobaric 
surface, the maximum permitted luminosity of a star is reduced by rotation, with 
respect to the usual Eddington limit. In the above relation, is the largest 
value of the opacity on the surface of the rotating star. For O-type stars with 
photospheric opacities dominated by electron scattering, the opacity k is the 
same everywhere on the star. 

For zero rotation, the usual expressions are found: = T — 47r *g M and 

t 4-rrcGM 
-^Edd = ~ ■ 

Critical limits. 

Critical limits correspond to values of respectively the luminosity and/or the 
velocity which impose that the total gravity become equal to zero at least at 
some peculiar places at the surface. We may identify different limits [39] : 

- We speak of the Eddington or P-limit, when rotation effects can be neglected 
and g r ad + 9grav = dE which implies that 
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<7grav is the gravitational acceleration. 
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In that case L = LEdd = ^ttcGM j ' n. The opacity k considered here is the total 
opacity. 

-The critical velocity or J7~limit is reached for a star with an angular velocity Q 
at the surface, when the effective gravity <j e ff = g sra _ v +g ro t = and in addition 
when radiation pressure effects can be neglected. 

-The fir-limit is reached when the total gravity gtot = 0, with significant effects 
of both rotation and radiation. This is the general case. It should lead to the two 
above cases in their respective limits. 

Using relation [Til we mav write the expression for the total gravity as 

<7tot = 9eff [1 - . (14) 

This shows that the expression for the total acceleration in a rotating star is sim- 
ilar to the usual one, except that T is replaced by the local value /^(t?). Indeed, 
contrarily to expressions such as gtot = <7eff (1 — r) often found in literature, 
we see that the appropriate Eddington factor given by Eq. [TTJ also depends on 
the angular velocity Q on the isobaric surface. 

From Eq. [TTJ we note that over the surface of a rotating star, which has a 
varying gravity and T e ff, is the highest at the latitude where is the 

largest, (if we neglect the effects of which is justified in general). If the 

opacity increases with decreasing T as in hot stars, the opacity is the highest at 
the equator and there the limit rh(i?) = 1 may by reached first. Thus, it is to 
be stressed that if the limit Tnld) — 1 happens to be met at the equator, it is 
not because g e ff is the lowest there, but because the opacity is the highest ! The 
reason for no direct dependence on £j e ff is because both terms g e ff cancel each 
other in the expression [TT1 of the flux ratio. 

The critical velocity is reached when somewhere on the star one has gtot = 0, 

i.e. 



fleff [l-r O (0)]=O. (15) 

This equation has two roots. The first one v cr it,i is given by the usual condition 
9eff = 0, which implies the equality fl 2 R^ h /(GM) = 1 at the equator. This 
corresponds to an equatorial critical velocity 

/ 2 GM \ 5 , . 

«crit,l = " -Rcb = ^3"^ j ■ ( 16 ) 

i? e b and R p b are respectively the equatorial and polar radius at the critical ve- 
locity. We notice that the critical velocity v C rit,i is independent on the Eddington 
factor. To this extent, this is in agreement with [17j . The basic physical reason 
for this independence is quite clear: the radiative flux tends toward zero when 
the effective gravity is zero, thus there is no effect of the radiative acceleration 
when this occurs. 

Equation 1151 has a second root, which is given by the condition i~h(#) = 1. 
If we call the Eddington ratio -T ma x the maximum value of k($)L(P) / (AircGM) 
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over the surface (in general at equator), has shown that for r meiK < 0.639, 
no value of fl can lead to Fq{$) — 1. Physically this means that when the star is 
sufficiently far from the Eddington limit, the effects of rotation on the radiative 
equilibrium are not sufficient for lowering the Eddington luminosity such that 
it may have an impact on the value of the critical velocity. In that case, Eq. [15] 
has only one root given by the classical expression. 

For a given large enough -T max (i.e. larger than 0.639), a second root is ob- 
tained given by 



9 a 1-T RI(uj) 



- v, 



4 " criM V'(w) R? 



(17) 



P b 



The quantity to is the fraction J?/J? c of the angular velocity at break-up. 
The quantity V'(u) is the ratio of the actual volume of a star with rotation u> 
to the volume of a sphere of radius i? p b- V(u>) is obtained by the integration of 
the solutions of the surface equation for a given value of the parameter u. 

For r max > 0.639, this second root is inferior to the first one. Thus it is 
encountered first and is therefore the expression of the critical velocity that has 
to be used. 

The mass loss rates 

Due to the von-Zeipel theorem, the radiative flux, which is the driving force 
for the stellar winds of massive hot stars, varies as a function of the colatitude. 
This effect, when accounted for in the theory of the line driven wind theory, lead 
to an enhancement of the quantities of mass lost and to wind anisotropics. We 
shall describe in more details these effects in the Section 1.3 below. 

Three remarks 

1) The classical critical angular velocity or the J7-limit (to distinguish it from 
the J7_T-limit as defined by 39J) in the frame of the Roche model is given by 



where i? p b is the polar radius when the surface rotates with the critical velocity. 
The critical velocity is given by 



2 GM\ 

3 R p b J 



(19) 



Replacing fl in Eq.|4]by (J?/J2 C rit) ■ ^crit and using Eq.HH one obtains a relation 
between i? p , / and i? p b, 

3 /^crit\ 2/3 /2(/ - 1)\ 1/3 , , 
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Thus one has 

With a good approximation (see below) one has that i? p b /Rp — 1 and therefore 

This equation is quite useful since it allows the determination of J? C rit from 
quantities obtained with a model computed for Q. This is not the case of Eq. [20] 
which involves R p h whose knowledge can only be obtained by computing models 
at the critical limit. In general Eq. 1221 gives a very good approximation of J? C rit 
(see [16] for a discussion of this point). 

Setting v the velocity at the equator, one has that 

v _ QR C _ H R c Rp R p b ^3-) 



Wcrit ^cRch ^crit Rp Rpb Reb 

where R c b is the equatorial radius when the surface rotates with the critical 
velocity. Using Eq. [20] above, and the fact that in the Roche model R p b/R e b = 
2/3, one obtains 

n ^ 1/3 



^crit V 



2(/ - 1) (24) 



The relations between u/w C rit and J?/]7 cl -it obtained in the frame of the Roche 
model (see Eq. |2~4"]) for the 1 and 60 M Q stellar models at Z — 0.02 are shown in 
Fig.[T] In case we suppose that the polar radius remains constant i? p b/i? p = 1, 
then Eq. [221 can be used and one obtains a unique relation between J7/i7 cr i t and 
u/fcrit, independent of the mass, metallicity and evolutionary stage considered. 
One sees that the values of u/fcrit is smaller than that of Q/ Q cv n by at most 
~25%. At the two extremes the ratios are of course equal. 

An interesting quantity is the ratio of the centrifugal acceleration, a ccn , to 
the gravity, g , at the equator 



Aeon _ ^ 2 i?c _ / & \ 2 /2\ 3 .3 / -Rp 



g c GM \O CTit J \3J f \Rp h ) ' (25) 

where we have used Eq. [TBI and divided/multiplied by Rp h . Replacing R p b/R p 
by its expression deduced from Eq. 1201 we obtain 

^=2(/-l). (26) 
9c 

We can check that at the critical limit, when / = 3/2, then a con = 9c- 

2) It is interesting to note that the fact that the effective temperature varies 
as a function of the colatitude on an isobaric surface does not necessarily imply 
that the temperature varies as a function of the colatitude on an isobaric surface. 
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Fig. 1. Relation between v/v CT i t and Q/ Q cr it obtained in the frame of the Roche model. 
The continuous line is obtained assuming i? p b/-R P = 1 (see text and Eqs. 1221 and I24|! . 
The dot-dashed lines show the relations for the Z=0.02 models with 1 and 60 Mq using 
Eg. 1211 The dotted line is the line of slope 1. 

For instance in the conservative case, the temperature is constant on equipo- 
tentials which are also isobaric surfaces. This simply illustrates the difference 
between the effective temperature whose definition is related to the radiative 
flux (F = crT^ff) and hence to the temperature gradient and the temperature 
itself. Interestingly, one has that in a conservative case, Fq{&) is constant on iso- 
baric surfaces (neglecting the term CC$))- This means that when the J?.T-limit 
is reached, it is reached over the whole stellar surface at the same time. This is 
in contrast with the f?-limit which is reached first at the equator. 

3) Recently we have examined the effects of rotation on the thermal gradient 
and on the Solberg-Hoiland term by analytical developments and by numerical 
models [35]. Writing the criterion for convection in rotating envelopes, we show 
that the effects of rotation on the thermal gradient are much larger and of op- 
posite sign to the effect of the Solberg-Hoiland criterion. On the whole, rotation 
favors convection in stellar envelopes at the equator and to a smaller extent 
at the poles. In a rotating 20 M star at 94% of the critical angular velocity, 
there are two convective envelopes, the biggest one has a thickness of 13.2% of 
the equatorial radius. The convective layers are shown in Fig. [5] They are more 
extended than without rotation. In the non-rotating model, the corresponding 
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Fig. 2. 2-D representation of the external convective zones (in red) and of the con- 
vective core (blue) in a model of 20 Mq with X = 0.70 and Z — 0.020 at the end of 
MS evolution with fast rotation ( f2/O CI it = 0.94). The axis are in units of cm. Figure 
taken from [38] . 

convective zone has a thickness of only 4.6% of the radius. The occurrence of 
outer convection in massive stars has many consequences (see |38j). 

1.2 Transport mechanisms of angular momentum and of chemical 
species 

In a solid body rotation state, the first instability to set up is a thermal instabil- 
ity called the meridional circulation (see below). It consists in large meridional 
currents which transport angular momentum either inside-out or conversely from 
the outer regions toward the inner ones. Such meridional currents rapidly build 
ut gradients of the angular momentum both in the "horizontal direction" (i.e. 
along isobaric surface) and in the vertical one. Along isobaric surface, any gra- 
dient of fl triggers a strong horizontal turbulence. Indeed in that direction the 
instability can develop without having to overcome any stable density gradient. 
As a consequence any gradient of fl along isobaric surfaces is rapidly erased and 
the star settles into a "shellular" rotation state [64] . This means that fl can be 
considered as nearly constant along isobars. In the vertical direction, where in 
a radiative zone, a stable density gradient counteracts any instability, the gra- 
dients of fl are eroded on much longer timescales (see below). The equations 
below describe the interactions of meridional currents and of shear instabilities 
in a state of shellular rotation [64] . 
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Meridional circulation. Meridional circulation is an essential mixing mecha- 
nism in rotating stars and there is a considerable literature on the subject (see 
ref. in [57]). The velocity of the meridional circulation in the case of shellular 
rotation was derived by [64] . The velocity of meridional circulation is derived 
from the equation of energy conservation [43] 



P T 



dS , . 

Tt + {err 



U) • vs 



= div(xVT) + pe - divF,, 



(27) 



where S is the entropy per unit mass, \ the thermal conductivity, e the rate of 
nuclear energy per unit mass and the flux of thermal energy due to horizontal 
turbulence. All the quantities are expanded linearly around their average on a 
level surface or isobar, using Legendre Polynomials P2(cos9). For instance 

T{P,0)=T{P)+TP 2 (cos 6) . 

Then Eq. [27] is linearized and an expression for C/ 2 can be deduced [64] • Using 
the same method [36] revised the expression for U 2 to account for expansion 
and contraction in non-stationary models. They also studied the effects of the 
/^-gradients (mean molecular weight gradients) , of the horizontal turbulence and 
considered a general equation of state. They obtained 



U 2 (r) 



P 



pgCpTlV^ - V + (v>/<$)V M ] 



—{En 



CpdO 
~T~dt 



(28) 



where M± = Mil 



is the reduced mass and the other symbols have 
and [36J 0. The driving term in the square brackets 



1 



Q 



2ttGp 



( n 2 r 3 \ 

\ GM J 



2ttGp„ 

the same meaning as in [BJ 

in the second member is Eq ■ It behaves mainly like Eq 

The term p means the average on the considered equipotential. The term with 
the minus sign in the square bracket is the Gratton-Opik term, which becomes 
important in the outer layers when the local density is small. This term produces 
negative values of Ua(r) (noted U(r) from now), meaning that the circulation is 
going down along the polar axis and up in the equatorial plane. This makes an 
outward transport of angular momentum, while a positive U(r) gives an inward 
transport. At lower Z, the Gratton-Opik term is negligible, which contributes 
to make larger /^-gradients in lower Z stars. 

Recently [41J rederived the system of partial differential equations, which 
govern the transport of angular momentum, heat and chemical elements. They 
expand the departure from spherical symmetry to higher order and include ex- 
plicitly the differential rotation in latitude, to first order. Boundary conditions 
for the surface and at the frontiers between radiative and convective zones are 
also explicitly given in this paper. 



Shellular rotation. The differential rotation which results from the evolution 
and transport of the angular momentum makes the stellar interior highly tur- 
bulent. As explained above, the turbulence is very anisotropic, with a much 

3 U2 is the same as U in Eq. [51] 
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stronger geostrophic-like transport in the horizontal direction than in the verti- 
cal one [64], where stabilisation is favoured by the stable density gradient. This 
strong horizontal transport is characterized by a large diffusion coefficient D^. 
Various expressions have been proposed: 

• The usual expression for the coefficient of viscosity due to horizontal 
turbulence and for the coefficient D^ of horizontal diffusion, which is of the 
same order, is, according to [64] . 

D h ~ vy, = —r \2V(r) - aU(r)\ , (29) 

Ch 

where r is the appropriately defined eulerian coordinate of the isobar [44 . 
V(r) is defined by ug(r, 6) = V(r) dP2< ^° s9 ^ where uo is the horizontal com- 
ponent of the velocity of the meridional current^, a = h d ^ — and Ch is 
a constant of order of unity or smaller. This equation was derived assuming 
that the differential rotation on an isobaric surface is small [36 . 

• [29] has derived an expression for the coefficient of diffusion by horizontal 
turbulence in rotating stars. He has obtained 

D h (x r (rY3(r)V [2V-aU])3 . (30) 

This expression can be written in the usual form — | I • v for a viscosity, 
where the appropriate velocity v is a geometric mean of 3 relevant velocities: 
a velocity (2V — all) as in Eq[29]by [64], the horizontal component V of 
the meridional circulation, the average local rotational velocity rQ(r). This 
rotational velocity is usually much larger than either U(r) or V(r), typically 
by 6 to 8 orders of a magnitude in an upper Main Sequence star rotating 
with the average velocity. 

• From torque measurements in the classical Couette-Tayler experiment [52] , 
[42] have found the following expression 

" h= {lo) 1/2 H2V-aC/|])^ , (31) 

with « 1.5 x 1CT 5 [52]. 

The horizontal turbulent coupling favours an essentially constant angular veloc- 
ity J? on the isobars. This rotation law, constant on shells, applies to fast as well 
as to slow rotators. As an approximation, it is often represented by a law of the 
form Q = fi(r) ([64]: see also [I5J). Let us note here that the exact value of the 
diffusion coefficient Dh is not well known. Indeed the values of the numerical 
factors intervening in the various expressions shown above may vary to some 
extent. Since the expression of Dh intervenes in the formulas for U r , for D s hear 
and D c g (see below), these uncertainties have some impact on the amplitudes 
of the transport mechanisms. 



4 V(r) can be obtained from U(r) see Eq. 2.10 in [64] 
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Shear turbulence and mixing. In a radiative zone, shear due to differential 
rotation is likely to be a most efficient mixing process. Indeed shear instability 
grows on a dynamical timescale that is of the order of the rotation period [64] . 
The usual criterion for shear instability is the Richardson criterion, which com- 
pares the balance between the restoring force of the density gradient and the 
excess energy present in the differentially rotating layers, 

A 2 , 1 
Ri = ad .. n < -, (32) 

where we have taken the average over an isobar, r is the radius and A a( j the 
Brunt- Vaisala frequency given by 



gs_ 

Hp 



ad 



(33) 



When thermal dissipation is significant, the restoring force of buoyancy is re- 
duced and the instability occurs more easily. Its timescale is however longer, 
being the thermal timescale. This case is referred to as "secular shear instabil- 
ity". The criterion for low Peclet numbers Pe (i.e. of large thermal dissipation, 
see below) has been considered by [63] , while the cases of general Peclet numbers 
Pe have been considered by [27], [32], who give 



Ri = 



(0.8836 f2 



dlnj]_\2 
d In r 



r 



r + i 



(v ad - V) 



< 



(34) 



The quantity r = Pe/6, where the Peclet number Pe is the ratio of the thermal 
cooling time to the dynamical time, i.e. Pe — where v and £ are the char- 
acteristic velocity and length scales, and K = (4acT 3 ) / (3CpKp 2 ) is the thermal 
diffusivity. A discussion of shear-driven turbulence by [5] suggests that the lim- 
iting Ri number may be larger than 4. 

To account for shear transport and diffusion, we need a diffusion coefficient. 
Amazingly, a great variety of coefficients D s hear 
applied (see a more extended discussion in [45]): 



^v£ have been derived and 



1. _64j defines the diffusion coefficient corresponding to the eddies which have 
the largest Pe number so that the Richardson criterion is just marginally 
satisfied. However, the effects of the vertical /x-gradient are not accounted for 
and the expression only applies to low Peclet numbers. The same has been 
done by [52) . who considers also the effect of the vertical ^(-gradient, the 
case of general Peclet numbers and, in addition they account for the coupling 
due to the fact that the shear also modifies the local thermal gradient. This 
coefficient has been used by [44] and by [T5J. The comparisons of model 
results and observations of surface abundances have led many authors to 
conclude that the /x-gradients appear to inhibit the shear mixing too much 
with respect to what is required by the observations ( [7] . [44] . [T9] ) . 
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2. Instead of using a gradient V M in the criterion for shear mixing, [7] and 
[19] write /^V^ with a factor / M = 0.05 or even smaller. This procedure is 
not satisfactory since it only accounts for a small fraction of the existing 
/^-gradients in stars. The problem is that the models depend at least as 
much (if not more) on / M than on rotation, i. e. a change of in the allowed 
range (between and 1) produces as important effects as a change of the 
initial rotational velocity. This situation has led to two other more physical 
approaches discussed below. Also [19] introduces another factor f c to adjust 
the ratio of the transport of the angular momentum and of the chemical 
elements like [49]. 

3. Around the convective core in the region where the /i-gradient inhibits mix- 
ing, there is anyway some turbulence due to both the horizontal turbulence 
and to the semiconvective instability, which is generally present in massive 
stars. This situation has led to the hypothesis [28] that the excess energy in 
the shear, or a fraction a of it of the order of unity, is degraded by turbulence 
on the local thermal timescale. This progressively changes the entropy gra- 
dient and consequently the /i-gradient. This hypothesis leads to a diffusion 
coefficient -Dshcar given by 



Nf 



K 

ad 



\a ( 0.8836 ) -(V'-V) 

4 V d\nr ' 



(35) 



The term V'— V in Eq.l3"5lexpresses either the stabilizing effect of the thermal 
gradients in radiative zones or its destabilizing effect in semiconvective zones 
(if any). When the shear is negligible, -D s hcar tends toward the diffusion 
coefficient for semiconvection by |25] in semiconvective zones. When the 
thermal losses are large (V = V), it tends toward the value 

(d In f2 \ ^ 
0.8836 fl-^) , (36) 

given by [64) . Eq.[35]is completed by the three following equations expressing 
the thermal effects [35] 

n ok r v Vrad + ( ^T )Vad (vf \ 
V she ^ = 2Kl V = 6r2 , (37) 

1 + li+rJ 



V'-V= -^(Vad-V). (38) 

The system of 4 equations given by Eqs.l3l^l3"7landl3"51form a coupled system 
with 4 unknown quantities Z? s hcar, r, V and V. The system is of the third 
degree in r. When it is solved numerically, we find that as a matter of fact 
the thermal losses in the shears are rather large in massive stars and thus 
that the Peclet number Pe is very small (of the order of 10~ 3 to 10 -4 ). For 
very low Peclet number Pe — 6T, the differences (V — V) are also very small 
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as shown by Eq. [351 Thus, we conclude that Eq. [35] is essentially equivalent, 
at least in massive stars, to the original Eq. [36] above, as given by [64]. We 
may suspect that this is not necessarily true in low and intermediate mass 
stars since there the Pe number may be larger. 
4. [55] found that the diffusion coefficient for the shears is modified by the 
horizontal turbulence. The change can be an increase or a decrease of the 
diffusion coefficient depending on the various parameters, as discussed below. 
Thus, we have 



where is the coefficient of horizontal diffusion (cf. [M])- We ignore here 
the thermal coupling effects discussed by Maeder ([28]) because they were 
found to be relatively small and they increase the numerical complexity. 
Interestingly, we see that in regions where V M ~ 0, Eq. [39] leads us to replace 
K by (K + Dh) in the usual expression (cf. [55]), i.e. it reinforces slightly 
the diffusion in regions which are close to chemical homogeneity. On the 
contrary, in regions where V M dominates with respect to (V a( j — V ra d), the 
transport is proportional to Dh rather than to K, which is quite logical 
since the diffusion is then determined by Z\ rather than by thermal effects. 
The above result shows the importance of the treatment for the meridional 
circulation, since in turn it determines the size of Z?h and to some extent the 
diffusion by shears. 

Of course, the Reynolds condition D s hoar > ^iyRe c must be satisfied in order 
that the medium is turbulent. The quantity v is the total viscosity (radiative + 
molecular) and Re c the critical Reynolds number estimated to be around 10 (cf. 
[13] : [64]). The numerical results indicate that the conditions for the occurrence 
of turbulence are satisfied. 

Transport of the angular momentum Let us express the rate of change 
of the angular momentum, of the clement of mass in the volume ABCD 
represented in Fig. [3] 



where M is the momentum of the forces acting on the volume element. We 
assume that angular momentum is transported only through advection (by a 
velocity field U) and through turbulent diffusion, which may be different in the 
radial (vertical) and tangential (horizontal) direction. The component of the 



D = 



(K + A,) 



(39) 



Hp 
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angular momentum aligned with the rotational axis is equal to 
pr 2 sin 9d9dipdr r sin 9Q y rsrnO, 

Mass of ABCD velocity distance to axis 

where fl = (p. Since the mass of the volume element ABCD does not change, 
the rate of change of the angular momentum can be written 

pr 2 sin 6d9dtpdr-^(r 2 sin 2 6Q) Mr . (40) 

Due to shear, forces apply on the surfaces of the volume element. The force on 
the surface AB is equal to 

dn 2 

Tj v rsmO-^— r sin^d^dt^. 

vertical viscosity v surface AB 

vertical shear 

The component of the momentum of this force along the rotational axis is 
r] v r sin 6-^-d6dtp r sin 0. 

v ■ distance to axis 

force on AB 

The component along the rotational axis of the resultant momentum of the forces 
acting on AB and CD is equal to 

8 3Q 

— (r] v r 4 sm 3 dd9dip—)dr. (41) 
or or 

The force on the surface AC due to the tangential shear is equal to 

r\h f sin 9 r sin 9d(pdr , 

v surface AC 

tangential shear 

where r]h is the horizontal viscosity. The component along the rotational axis of 
the resultant momentum of the forces acting on AC and BD is equal to 

— (rj h r 2 sin 3 edrd^p— )rd0. (42) 

Using Eqs. 201 ED an d HH simplifying by drd^dt^, one obtains the equation for 
the transport of the angular momentum 

pr 2 sin^ (r2 sin2 6 Q) Mr = l( n y sin 3 df) + ^r 2 sin 3 9^). (43) 



5 The components perpendicular to the rotational axis cancel each other when the 
integration is performed over ip. 
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Fig. 3. The momentum of the viscosity forces acting on the element ABCD is derived 
in the text and the general form of the equation describing the change with time of the 
angular momentum of this element is deduced. The star rotates around the vertical 
axis with the angular velocity Q\ r and 9 are the radial and colatitude coordinates of 
point A. 



Setting r] v = pD v and rjh — pDh and dividing the left and right member by 
r 2 sin#, one obtains 

d / 2 • 2 /i ^\ sm2 # ® , 4 <9i? , 1 d . _ . -> . dfl . , A1 . 

P -y W en )Mr = -^-ony-) + — elfe {pD h sm" e-). (44) 



Now, the left-handside term can be written 



p A( r 2 sin 2 eQ)Mr = iL (/9r 2 sin 2 6Q)m ^ _ r 2 sin 2 en^\ Mr . 



d_ 

dV 



Using the relation between the Lagrangian and Eulerian derivatives, one has 

P ^-{r 2 sin 2 en) Mr = 

at 



^-(pr 2 sin 2 en) r + U • V(pr 2 sin 2 00) - r 2 sin 2 9tt^- 
ot ' at 



M r - 



Using 

and the continuity equation 



dp dp 
dt^=Q- t \r + V.Vp, 



dp 

— \ r = -div(pU), 



one obtains dpjdt\M r + pdivU = 0, which incorporated in Eq. l45l gives 
p^(r 2 sin 2 6Q) Mr = ^(pr 2 sin 2 BQ) r + V(Upr 2 sin 2 9Q). 



(45) 
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Developing the divergence in spherical coordinates and using Eq. 1441 one finally 
obtains the equation describing the transport of the angular momentum ( [36] ; 

d 1 d 1 d 

— {pr 2 sin 2 9Q) r + ^^(pr 4 sin 2 9w r Q) + —^^{pr 2 sin 3 6w e Q) = 
at r z or r sin oO 



sm 



2 , 



d A dfl 1 d . 3 df2 



dr^ r ^ + ^ede^ Dh ^ e m^ (46) 

where w r = U r + r is the sum of the radial component of the meridional circu- 
lation velocity and the velocity of expansion/contraction, and wg = Ug, where 
Ug is the horizontal component of the meridional circulation velocity. Assuming, 
as in [M], that the rotation depends little on latitude due to strong horizontal 
diffusion, we write 

f2(r,d) = Q{r) + Q{r,6), 

with fl <C fl. The horizontal average /? is defined as being the angular velocity 
of a shell rotating like a solid body and having the same angular momentum as 
the considered actual shell. Thus 

^_ / Q sin 3 0d6 
~ Jsin 3 6»d6> 

Any vector field whose Laplacian is nul can be decomposed in spherical harmon- 
ics. Thus, the meridional circulation velocity can be written [41] 

U = £ Ui (r)fl (cos 6) e r + £ Vt (r) e e , 



where e r and eg are unit vectors along the radial and colatitude directions re- 
spectively. Multiplying Eq. [46jby sin#d6> and integrating it over 9 from to n, 
one obtains [36] 

l- t ^o> r . £±<smv - «b + ~ ■ < 47 » 

It is interesting to note that only the I = 2 component of the circulation is 
able to advect a net amount of angular momentum. As explained in [54j the 
higher order components do not contribute to the vertical transport of angular 
momentum. Note also that the change in radius r of the given mass shell is 
included in Eq. I47i which is the Eulerian formulation of the angular momentum 
transport equation. In its Lagrangian formulation, the variable r is linked to 
M r through dM r = 4irr 2 pdr, and the equation for the transport of the angular 
momentum can be written 
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The characteristic time associated to the transport of f2 by the circulation is 

m 

ta^t K H\—-\ , (49) 

where g s is the gravity at the surface and t^H the Kelvin-Helmholtz timescale, 
which is the characteristic timescale for the change of r in hydrostatic models. 
From Eq. [49j one sees that tn is a few times txH , which itself is much shorter 
that the Main Sequence lifetime. 

For shellular rotation, the equation of transport of angular momentum in the 
vertical direction is in lagrangian coordinates (cf. [64]; [36] ) 



Q(r) is the mean angular velocity at level r. The vertical component u(r, 9) of 
the velocity of the meridional circulation at a distance r to the center and at a 
colatitudc 9 can be written 



u(r,9) = U(r)P 2 (cos 9), (51) 

where P2(cos9) is the second Legendre polynomial. Only the radial term U(r) 
appears in Eq. 1501 The quantity D is the total diffusion coefficient representing 
the various instabilities considered and which transport the angular momentum, 
namely convection, semiconvection and shear turbulence. As a matter of fact, 
a very large diffusion coefficient as in convective regions implies a rotation law 
which is not far from solid body rotation. In this work, we take D = Z? s hear in 
radiative zones, since as extra-convective mixing we consider shear mixing and 
meridional circulation. 

In case the outward transport of the angular momentum by the shear is com- 
pensated by an inward transport due to the meridional circulation, we obtain the 
local conservation of the angular momentum. We call this solution the stationary 
solution. In this case, U(r) is given by (cf. [64] ) 

The full solution of Eq. taking into account U(r) and D gives the non- 
stationary solution of the problem. In this case, fi(r) evolves as a result of the 
various transport processes, according to their appropriate timescales, and in 
turn differential rotation influences the various above processes. This produces 
a feedback and, thus, a self-consistent solution for the evolution of fl(r) has to 
be found. 




Fig.|3]shows the evolution of U(r) in a model of a 20 M star with Z — 0.004 
and an initial rotation velocity Wi n ; = 300 km s~ x [33]. U(r) is initially positive 
in the interior, but progressively the fraction of the star where U(r) is negative 
is growing. This is due to the Gratton-Opik term in Eq. (2), which favors a 
negative U(r) in the outer layers, when the density decreases. This negative 
velocity causes an outward transport of the angular momentum, as well as the 
shearqS 

6 When U is negative, the meridional currents turn anticlockwise, i.e. go inwards along 
directions parallel to the rotational axis and go outwards in directions parallel to the 
equatorial plane. 
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The transport of angular momentum by circulation has often been treated 
as a diffusion process ([15]: [49] : [19]). From Eq. [50] we see that the term with U 
(advection) is functionally not the same as the term with D (diffusion). Physi- 
cally advection and diffusion are quite different: diffusion brings a quantity from 
where there is a lot to other places where there is little. This is not necessarily 
the case for advection. A circulation with a positive value of U(r), i.e. rising 
along the polar axis and descending at the equator, is as a matter of fact making 
an inward transport of angular momentum. Thus, we see that when this process 
is treated as a diffusion, like a function of even the sign of the effect may be 
wrong. 

The expression of U(r) given above (Eq. [25]) involves derivatives up to the 
third order, thus Eq. [50] is of the fourth order, which makes the system very 
difficult to solve numerically. In practice, we have applied a Henyey scheme to 
make the calculations. Eq. [50] also implies four boundary conditions. At the 
stellar surface, we take (cf. [56] ) 



— = and U(r) = (53) 



and at the edge of the core we have 



— = and Q(r) = Q COIC . (54) 

We assume that the mass lost by stellar winds is just embarking its own angular 
momentum. This means that we ignore any possible magnetic coupling, as it 
occurs in low mass stars. It is interesting to mention here, that in case of no 
viscous, nor magnetic coupling at the stellar surface, i.e. with the boundary 
conditions [53] the integration of Eq. [50] gives for an external shell of mass AM 



AM^(f2r 2 ) = -—pr 4 f2U{r). (55) 
dt 5 

This equation is valid provided the stellar winds are spherically symmetric. When 
the surface velocity approches the critical velocity, it is likely that there are 
anisotropies of the mass loss rates (polar ejection or formation of an equato- 
rial ring) and thus the surface condition should be modified according to the 
prescriptions of [25] . 



Mixing and transport of the chemical elements. A diffusion-advection 
equation like Eq. l50l should normally be used to express the transport of chemical 
elements. However, if the horizontal component of the turbulent diffusion is 
large, the vertical advection of the elements can be treated as a simple diffusion 
[6] with a diffusion coefficient -D c ff- As emphasized by [Bj , this does not apply 
to the transport of the angular momentum. D e g is given by 
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D CS = , (56) 

where is the coefficient of horizontal turbulence. Eq. [56] expresses that the 
vertical advection of chemical elements is severely inhibited by the strong hor- 
izontal turbulence characterized by D^. Thus, the change of the mass fraction 
Xi of the chemical species i is simply 



The second term on the right accounts for composition changes due to nuclear 
reactions. The coefficient D mlK is the sum D mlx = D s h ea r + Defi and D e ff is 
given by Eq. 1561 The characteristic time for the mixing of chemical elements 
is therefore i m i x ~ -jp— and is not given by t c - ltc ~ ^, as has been generally 
considered |53j . This makes the mixing of the chemical elements much slower, 
since D c g is very much reduced. In this context, we recall that several authors 
have reduced by large factors, up to 30 or 100, the coefficient for the transport of 
the chemical elements, with respect to the transport of the angular momentum, 
in order to better fit the observed surface compositions (cf. [IS]). This reduction 
of the diffusion of the chemical elements is no longer necessary with the more 
appropriate expression of D c g given here. 

Surface enrichments due to rotation are illustrated in Fig [5] The tracks are 
plotted in the plane (N/C)/(N/C); n i versus P where P is the rotational period in 
hours. During the evolution the surface is progressively enriched in CNO burning 
products, i.e. is enriched in nitrogen and depleted in carbon. At the same time, 
the rotational period increases. 

When the effects of the shear and of the meridional circulation compensate 
each other for the transport of the angular momentum {stationary solution), the 
value of U entering the expression for D e g is given by Eq. [52l 



1.3 Rotation and mass loss 

We can classify the effects of rotation on mass loss in three categories. 

1. The structural effects of rotation. 

2. The changes brought by rotation on the radiation driven stellar winds. 

3. The mass loss induced by rotation at the critical limit. 

Let us now consider in turn these various processes. 



Structural effects of rotation on mass loss. Rotation, by changing the 
chemical structure of the star, modifies its evolution. For instance, moderate 
rotation at metallicities of the Small Magellanic Cloud (SMC) favors redward 
evolution in the Hertzsprung-Russel diagram. This behavior can account for the 
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Fig. 5. Evolutionary tracks in the plane surface N/C ratio, normalized to its initial 
value, versus the rotational period in hours for different initial mass stars, various initial 
velocities and for the metallicities Z=0.02 and 0.002. Positions of some periods in days 
are indicated at the bottom of the figure. The dotted tracks never reach the critical 
limit during the MS phase. The short dashed tracks reach the critical limit during the 
MS phase. The dividing line between the shaded and non-shaded areas corresponds to 
the entrance into the phase when the star is at the critical limit during the MS phase. 
If Be stars are stars rotating at or very near the critical limit, present models would 
predict that they would lie in the vicinity of this dividing line or above it. Note the 
different vertical scales used when comparing similar masses at different metallicities. 
Figure taken from [16] . 
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high number of red supergiants observed in the SMC [33 , an observational fact 
which is not at all reproduced by non-rotating stellar models. 

Now it is well known that the mass loss rates are greater when the star evolves 
into the red part of the HR diagram, thus in this case, rotation modifies the mass 
loss indirectly, by changing the evolutionary tracks. The v in i = 0, 200, 300 and 
400 km s _1 models lose respectively 0.14, 1.40, 1.71 and 1.93 M© during the core 
He-burning phase (see Table 1 in [33]). The enhancement of the mass lost reflects 
the longer lifetimes of the red supergiant phase when velocity increases. Note 
that these numbers were obtained assuming that the same scaling law between 
mass loss and metallicity as in the MS phase applies during the red supergiant 
phase. If, during this phase, mass loss comes from continuum-opacity driven wind 
then the mass-loss rate will not depend on metallicity (see the review by [58]). 
In that case, the redward evolution favored by rotation would have a greater 
impact on mass loss than that shown by the computations shown above. 

Of course, such a trend cannot continue forever. For instance, at very high 
rotation, the star will have a homogeneous evolution and will never become a 
red supergiant [2B]. In this case, the mass loss will be reduced, although this 
effect will be somewhat compensated by other processes: first by the fact that 
the Main-Sequence lifetime will last longer, second, by the fact that the star 
will enter the Wolf-Rayet phase (a phase with high mass loss rates) at an earlier 
stage of its evolution, and third by the fact that the star may encounter the 
Q- limit. 

Radiation driven stellar winds with rotation. The effects of rotation on 
the radiation driven stellar winds result from the changes brought by rotation to 
the stellar surface. They induce changes of the morphologies of the stellar winds 
and increase their intensities. 

Stellar wind anisotropics 

Naively we would first guess that a rotating star would lose mass preferen- 
tially from the equator, where the effective gravity (gravity decreased by the ef- 
fect of the centrifugal force) is lower. This is probably true when the star reaches 
the f?-limit (i.e. when the equatorial surface velocity is such that the centrifugal 
acceleration exactly compensates the gravity), but this is not correct when the 
star is not at the critical limit. Indeed as recalled above, a rotating star has a non 
uniform surface brightness, and the polar regions are those which have the most 
powerful radiative flux. Thus one expects in case the opacity does not vary at 
the surface, that the star will lose mass preferentially along the rotational axis. 
This is correct for hot stars, for which the dominant source of opacity is electron 
scattering. In that case the opacity only depends on the mass fraction of hydro- 
gen and does not depends on other physical quantities such as temperature. In 
that way, rotation induces anisotropies of the winds ([31]; [14]). This is illustrated 
in the left panel of Fig. [6] Wind anisotropies have consequences for the angular 
momentum that a star retains in its interior. Indeed, when mass is lost preferen- 
tially along the polar axis, little angular momentum is lost. This process allows 
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loss of mass without too much loss of angular momentum a process which might 
be important in the context of the evolutionary scenarios leading to Gamma 
Ray Bursts. Indeed in the framework of the collapsar scenario ([62]), one has to 
accommodate two contradictory requirements: on one side, the progenitor needs 
to lose mass in order to have its H and He-rich envelope removed at the time 
of its explosion, and on the other hand it must have retained sufficient angular 
momentum in its central region to give birth to a fast rotating black-hole. 

Intensities of the stellar winds 

The quantity of mass lost through radiatively driven stellar winds is enhanced 
by rotation. This enhancement can occur through two channels: by reducing the 
effective gravity at the surface of the star, by increasing the opacity of the outer 
layers through surface metallicity enhancements due to rotational mixing. 

• reduction of the effective gravity: The ratio of the mass loss rate of a star 
with a surface angular velocity fl to that of a non-rotating star, of the same 
initial mass, metallicity and lying at the same position in the HR diagram 
is given by [3S] 

(1 - r) " X -, (58) 



M(0) 



9 V fcrft.i ' 



where F is the electron scattering opacity for a non-rotating star with the 

same mass and luminosity, a is a force multiplier [24] . The enhancement 

factor remains modest for stars with luminosity sufficiently far away from the 

M(n) 



Eddington limit |39j . Typically, ; ~ 1.5 for main-sequence B-stars. In 
that case, when the surface velocity approaches the critical limit, the effective 
gravity decreases and the radiative flux also decreases. Thus the matter 
becomes less bound when, at the same time, the radiative forces become 
also weaker. When the stellar luminosity approaches the Eddington limit, 
the mass loss increases can be much greater, reaching orders of magnitude. 
This comes from the fact that rotation lowers the maximum luminosity or 
the Eddington luminosity of a star. Thus it may happen that for a velocity 
still far from the classical critical limit, the rotationally decreased maximum 
luminosity becomes equal to the actual luminosity of the star. In that case, 
strong mass loss ensues and the star is said to have reached the flF limit 

Effects due to rotational mixing: During the core helium burning phase, at 
low metallicity, the surface may be strongly enriched in both H-burning and 
He-burning products, i.e. mainly in nitrogen, carbon and oxygen. Nitrogen 
is produced by transformation of the carbon and oxygen produced in the He- 
burning core and which have diffused by rotational mixing in the H-burning 
shell [46] . Part of the carbon and oxygen produced in the He-core also diffuses 
up to the surface. Thus at the surface, one obtains very high value of the 
CNO elements. For instance a 60 M Q with Z=10 -8 and Uj n i = 800 km s _1 
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has, at the end of its evolution, a CNO content at the surface equivalent to 1 
million times its initial metallicity! In case the usual scaling laws linking the 
surface metallicity to the mass loss rates is applied, such a the star would 
lose due to this process more than half of its initial mass. 

Mass loss induced by rotation As recalled above, during the Main-Sequence 
phase the core contracts and the envelope expands. In case of local conservation 
of the angular momentum, the core would thus spin faster and faster while the 
envelope would slow down. In that case, it can be easily shown that the surface 
velocity would evolve away from the critical velocity (see e.g. [47]). In models 
with shellular rotation however an important coupling between the core and the 
envelope is established through the action of the meridional currents. As a net 
result, angular momentum is brought from the inner regions to the outer ones. 
Thus, would the star lose no mass by radiation driven stellar winds (as is the 
case at low Z), one expects that the surface velocity would increase with time 
and would approach the critical limit. In contrast, when radiation driven stellar 
winds are important, the timescale for removing mass and angular momentum 
at the surface is shorter than the timescale for accelerating the outer layers by 
the above process and the surface velocity decreases as a function of time. It 
evolves away from the critical limit. Thus, an interesting situation occurs: when 
the star loses little mass by radiation driven stellar winds, it has more chance to 
lose mass by reaching the critical limit. On the other hand, when the star loses 
mass at a high rate by radiation driven mass loss then it has no chance to reach 
the critical limit and thus to undergo a mechanical mass loss. This is illustrated 
in the right panel of Fig. [6] 

Discussion. At this point it is interesting to discuss three aspects of the various 
effects described above. First what are the main uncertainties affecting them? 
Second, what are their relative importance? And finally what are their conse- 
quences for the interstellar medium enrichment? 

Uncertainties 

In addition to the usual uncertainties affecting the radiation driven mass loss 
rates, the above processes poses three additional problems: 

1. What does happen when the CNO content of the surface increases by six 
orders of magnitude as was obtained in the 60 Mq model described above? 
Can we apply the usual scaling law between Z and the mass losses? This 
is what we have done in our models, but of course this should be studied 
in more details by stellar winds models. For instance, for WR stars, [59] 
have shown that at Z = Z Q /30, 60% of the driving is due to CNO elements 
and only 10% to Fe. Here the high CNO surface enhancements result from 
rotational mixing which enrich the radiative outer region of the star in these 
elements, but also from the fact that the star evolves to the red part of the HR 
diagram, making an outer convective zone to appear. This convective zone 
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Fig. 6. Left panel: Iso-mass loss distribution for a 120 Mq star with Log L/L0=6.O 
and T e ff = 30000 K rotating at a fraction 0.8 of critical velocity (figure from 31). 
Right panel: Evolution of the surface velocities for a 60 Mq star with 3 different initial 
metallicities. 



plays an essential role in dredging up the CNO elements at the surface. Thus 
what is needed here is the effects on the stellar winds of CNO enhancements 
in a somewhat red part of the HR diagram (typical effective temperatures 
of the order of Log T g ~3.8). 
2. Do stars can reach the critical limit? For instance, [2] obtain that during pre- 
main sequence evolution of rapidly rotating massive stars, "equatorial mass 
loss" or "rotational mass ejection" never occur (see also [3]). In these mod- 
els the condition of zero effective gravity is never reached. However, these 
authors studied pre-main sequence evolution and made different hypothe- 
ses on the transport mechanisms than in the present work. Since they were 
interested in the radiative contraction phase, they correctly supposed that 
"the various instabilities and currents which transport angular momentum 
have characteristic times much longer than the radiative-contraction time" . 
This is no longer the case for the Main-Sequence phase. In our models, we 
consistently accounted for the transport of the angular momentum by the 
meridional currents and the shear instabilities. A detailed account of the 
transport mechanisms shows that they are never able to prevent the star 
from reaching the critical velocity. Another difference between the approach 
in the work of [3j and ours is that [2] consider another distribution of the 
angular velocity than in our models. They supposed constant Q on cylin- 
drical surface, while here we adopted, as imposed by the theory of [64] , a 
"shellular rotation law" . They resolved the Poisson equation for the gravi- 
tational potential, while here we adopted the Roche model. Let us note that 
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the Roche approximation appears justified in the present case, since only 
the outer layers, containing little mass, are approaching the critical limit. 
The majority of the stellar mass has a rotation rate much below the critical 
limit and is thus not strongly deformed by rotation. Thus these differences 
probably explain why in our models we reach situations where the effective 
gravity becomes zero. 
3. What does happen when the surface velocity reaches the critical limit? Let 
us first note that when the surface reaches the critical velocity, the energy 
which is still needed to make equatorial matter to escape from the potential 
well of the star is still important. This is because the gravity of the system 
continues of course to be effective all along the path from the surface to the 
infinity and needs to be overcome. If one estimates the escape velocity from 
the usual equation energy for a piece of material of mass m at the equator 
of a body of mass M, radius R and rotating at the critical velocity, 

1 2 I 2 GMm . , 

-mv CTit + -mv csc — = 0, (59) 

one obtains, using v^ lit = GM/R that the escape velocity is simply reduced 
by a factor l/v2 = 0.71 with respect to the escape velocity from a non- 
rotating body 0- Thus the reduction is rather limited and one can wonder if 
matter will be really lost. A way to overcome this difficulty is to consider the 
fact that, at the critical limit, the matter will be launched into a keplerian 
orbit around the star. Thus, probably, when the star reaches the critical 
limit an equatorial disk is formed like for instance around Be stars. Here we 
suppose that this disk will eventually dissipate by radiative effects and thus 
that the material will be lost by the star. 

Practically, in the present models, we remove the supercritical layers. This 
removal of material allows the outer layers to become again subcritical at 
least until secular evolution will bring again the surface near the critical limit 
(see [48] for more details in this process) . Secular evolution during the Main- 
Sequence phase triggers two counteracting effects: on one side, the stellar 
surface expands. Local conservation of the angular momentum makes the 
surface to slow down and the surface velocity to evolve away from the critical 
limit. On the other hand, meridional circulation continuously brings angular 
momentum to the surface and accelerates the outer layers. This last effect in 
general overcomes the first one and the star rapidly reach again the critical 
limit. How much mass is lost by this process? As seen above, the two above 
processes will maintain the star near the critical limit for most of the time. 
In the models, we adopt the mass loss rate required to maintain the star at 
about 95-98% of the critical limit. Such a mass loss rate is imposed as long as 
the secular evolution brings back the star near the critical limit. In general, 
during the Main-Sequence phase, once the critical limit is reached, the star 
remains near this limit for the rest of the Main-Sequence phase. At the end 
of the Main-Sequence phase, evolution speeds up and the local conservation 

7 We suppose here that the vector v CBC is normal to the direction of the vector u cr it. 
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of the angular momentum overcomes the effects due to meridional currents, 
the star evolves away from the critical limit and the imposed "critical" mass 
loss is turned off. 

2 "Spinstars" at very low metallicities? 

Let us call "spinstars" those stars with a sufficiently high initial rotation in 
order to have their evolution significantly affected by rotation. In this section, 
we present some arguments supporting the view according to which spinstars 
might have been more common in the first generations of stars in the Universe. 
A direct way to test this hypothesis would be to obtain measures of surface 
velocity of very metal poor massive stars and to see whether their rotation is 
superior to those measured at solar metallicity. At the moment, such measures 
can be performed only for a narrow range of metallicities for Z between 0.002 and 
0.020. Interestingly already some effects can be seen. For instance [22] presents 
measurements of the projected rotational velocities of a sample of 100 early B- 
type main-sequence stars in the Large Magellanic Cloud (LMC). He obtains that 
the stars of the LMC are more rapid rotators than their Galactic counterparts 
and that, in both galaxies, the cluster population exhibits significantly more 
rapid rotation than that seen in the field (a point also recently obtained by 
[2Tj). More recently [40] obtain that the angular velocities of B (and Be stars) 
are higher in the SMC than in the LMC and MW. 

For B-type stars, the higher values obtained at lower Z can be the result 
of two processes: 1) the process of star formation produces more rapid rotators 
at low metallicity; 2) the mass loss being weaker at low Z, less angular mo- 
mentum is removed from the surface and thus starting from the same initial 
velocity, the low Z star would be less slowed down by the winds. In the case 
of B-type stars, the mass loss rates are however quite modest and we incline to 
favor the first hypothesis, i.e. a greater fraction of fast rotators at birth at low 
metallicity. Another piece of argument supporting this view is the following: in 
case the mass loss rates are weak (which is the case on the MS phase for B-type 
stars), then the surface velocity is mainly determined by two processes, the ini- 
tial value on the ZAMS and the efficiency of the angular momentum transport 
from the core to the envelope. In case of very efficient transport, the surface will 
receive significant amount of angular momentum transported from the core to 
the envelope. The main mechanism responsible for the transport of the angular 
momentum is meridional circulation. The velocities of the meridional currents in 
the outer layers are smaller when the density is higher thus in more metal poor 
stars. Therefore, starting from the same initial velocity on the ZAMS, one would 
expect that B-type stars at solar metallicity (with weak mass loss) would have 
higher surface velocities than the corresponding stars at low Z. The opposite 
trend is observed. Thus, in order to account for the higher velocities of B-type 
stars in the SMC and LMC, in the frame of the present rotating stellar models, 
one has to suppose that stars on the ZAMS have higher velocities at low Z. Very 
interestingly, the fraction of Be stars (stars rotating near the critical velocity) 
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with respect to the total number of B stars is higher at low metallicity f[37j; 
[61]). This confirms the trend discussed above favoring a higher fraction of fast 
rotators at low Z. 

There are at least four other striking observational facts which might receive 
an explanation based on massive fast rotating models. First, indirect observa- 
tions indicate the presence of very helium-rich stars in the globular cluster wCen 
[50j. Stars with a mass fraction of helium, Y, equal to 0.4 seem to exist, together 
with a population of normal helium stars with Y = 0.25. Other globular clusters 
appear to host helium-rich stars [4] , thus the case of wCen is the most spectacu- 
lar but not the only one. There is no way for these very low mass stars to enrich 
their surface in such large amounts of helium and thus they must have formed 
from protostellar cloud having such a high amount of helium. Where does this 
helium come from? We proposed that it was shed away by the winds of metal 
poor fast rotating stars [341 . 

Second, in globular clusters, stars made of material only enriched in H- 
burning products have been observed (see the review by [18]). Probably these 
stars are also enriched in helium and thus this observation is related to the one 
reported just above. The difference is that proper abundance studies can be per- 
formed for carbon, nitrogen, oxygen, sodium, magnesium, lithium, fluorine . . . , 
while for helium only indirect inferences based on the photometry can be made. 
[11] propose that the matter from which the stars rich in H-burning products 
are formed, has been released by slow winds of fast rotating massive stars. Of 
course, part of the needed material can also be released by AGB stars. The 
massive star origin presents however some advantages: first a massive star can 
induce star formation in its surrounding, thus two effects, the enrichment and 
the star formation can be triggered by the same cause. Second, the massive star 
scenario allows to use a less flat IMF than the scenario invoking AGB stars [5T] . 
The slope of the IMF might be even a Salpeter's one in case the globular cluster 
lost a great part of its first generation stars by tidal stripping (see |12)V 

Third, the recent observations of the surface abundances of very metal poor 
halo stars^l show the need of a very efficient mechanism for the production of 
primary nitrogen [8]. As explained in [9]), a very nice way to explain this very 
efficient primary nitrogen production is to invoke fast rotating massive stars. 
Very interestingly, fast rotating massive stars help not only in explaining the 
behavior of the N/O ratio at low metallicity but also those of the C/O. Predic- 
tions for the behaviour of the 12 C/ 13 C ratios at the surface of very metal poor 
non-evolved stars have also been obtained [TO] , 

Fourth, below about [Fe/H] < -2.5, a significant fraction of very iron-poor 
stars are C-rich (see the review by [I]). Some of these stars show no evidence of 
s-process enrichments by AGB stars and are thus likely formed from the ejecta 
of massive stars. The problem is how to explain the very high abundances with 
respect to iron of CNO elements. [48] and [20] proposed that these stars might 
be formed from the winds of very metal poor fast rotating stars. It is likely that 

8 These stars are in the field and present [Fe/H] as low as -4, thus well below the 
metallicities of the globular clusters. 
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rotation also affects the composition of the ejecta of intermediate mass stars. [35] 
predict the chemical composition of the envelope of a 7 M Q E-AGB star which 
have been enriched by rotational mixing. The composition presents striking sim- 
ilarities with the abundance patterns observed at the surface of CRUMPS. The 
presence of overabundances of fluorine and of s-process elements might be used 
to discriminate between massive and intermediate mass stars. 

All the above observations seem to point toward the same direction, an im- 
portant population of spinstars at low Z. How many? What is the origin of the 
fast rotation? What are the consequences for the Gamma ray Burst progenitors? 
All these questions have still to be addressed in a quantitative way and offer nice 
perspective for future works. 
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